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2-D  Nonlinear  Theory  of  the  Free 
Electron  Laser  Amplifier  for  an 
Electron  Beam  with  Finite 
Axial  and  Transverse  Dimensions 

Many  current  experiments  of  the  free  electron  laser  (FEL),  ultilize  electron  beams  from  a 
millimeter  to  a  few  centimeters  in  pulse  length.  The  short  pulse  length  is  typical  of  high  energy 
accelerators  such  as  RF  Linacs  and  microtrons.  The  finite  length  effect  of  the  electron  beam  on 
the  radiation  was  found  to  be  important  in  the  Stanford  oscillator  experiment/0  Currently, 
many  experiments  in  the  amplifying  mode  are  being  conducted  with  the  short  electron  beam 
pulses.  We  have  an  analytical  expression  for  the  gain  pulse  of  the  radiation  field,  applicable  to 
these  experiments. 

One-dimensional  analysis  of  the  radiation  field  for  electron  beams  of  finite  length  have 
been  numerically  simulated  on  computers/2*40  The  effect  of  the  finite  transverse  dimensions 
was  either  not  included,  or  incorporated  through  filling  factors.  On  the  other  hand,  previous 
three-dimensions  self-consistent  formulation(5)  of  the  radiation  field  for  a  semi-infinitely  long 
electron  beam  in  the  amplifying  configuration  has  resulted  in  a  number  of  interesting  effects 
not  obtainable  by  the  1-D  formulation.  Numerical  effort  to  find  the  growth  of  the  3-D  radia¬ 
tion  field  on  the  finite  length  electron  beam  began  only  recently /6). 

In  this  paper,  we  will  present  a  fully  2-D,  self-consistent,  non-linear,  analytical  analysis  of 
the  FEL  process  in  the  amplifier  mode  of  operation  treating  the  finite  length  and  transverse 
effects  associated  with  both  the  electron  beam  and  the  radiation  beam.  Our  formulism  also 
includes  various  efficiency  enhancement  schemes:  (i)  contouring  in  the  longitudinal  direction 
the  amplitude  and/or  the  wavelength  of  the  magnetic  wiggler  field,  and  (ii)  applying  an  external 
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D.C.  electric  field.  Analytical  results  in  the  amplifying  configuration  are  obtained  in  the  low 
gain,  trapped  particle  regime. 

The  schematic  of  the  configuration  is  shown  in  Fig.  1.  The  generalized  vector  potentials 
of  the  right-handed,  helicd,  static  magnetic  wiggler  field  and  the  electromagnetic  radiation  field 
are 

A*  (z)  *  Ah.  (z)  [cos  (J*Q  kw  ( z ')  dz)  ex  +  sin  (  kw  (z')  dz')  ej  (\) 

\R  (jc,  y .  z,  /)  -  Ar  ( x .  y,  z,  t)  [cos  (  —  z  —  <*>  t  A  \p  (x  y,  z ,  t))  ex 

c 

-  sin  {—  z  -  a}  t  +  $  (x,  y,  z,  /))  ev]  (2) 

c 

where  Aw  and  kw  are  all  slowly  varying  amplitude  and  wave  number  of  the  wiggler  field  and  AR 
and  ^  are  slowly  varying  amplitude  and  phase  of  the  electromagnetic  radiation  field  following 
the  electron  pulse.  We  also  include  an  external  DC  electric  field,  £DC(z)  *  -d^Dc(z)/dz  ez 
for  the  purpose  of  efficiency  enhancement. 


In  this  analysis  we  will  not  consider  the  gradient  in  the  wiggler  field.  This  is  a  good 
approximation  if  kwrb  «  1,  where  rb  is  the  radius  of  the  electron  beam.  If  the  FEL  is  operat¬ 
ing  in  a  trapped  particle  mode,  we  also  require(5)  rb  <  (yz0  kwV]  (8  V2  yz(//8a),/2  ( Ar/Aw)va , 
where  yl0  ”  (1  -  vz20/c2)_l/2,  0a  -  kM»./(y0»»oC2),  y0  “  yroya,  ya  “ 

(1  +  |e|2  A2(0)/(mo  c4))l/2.  and  vz0  is  the  axial  velocity  at  z  ■>  0. 


The  electron  motion  can  be  described  in  terms  of  their  phase  4i  in  the  ponderomotive 
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where  i/i(x0,  y 0.  f0>  0  ”  J0  (k„(.z’)  +  ot/c)dz'  +  c ot  +  v>(x,  y,  ,  z.  t)  is  the  phase  for 


the  electron,  which  was  at  (jro.yK.fo)  at  t  -  0,  y  -  y,yL ,  y2  -  (1  -  v2/c2)~'n,  yL  - 
(1  +  |e|2/4»(r)/(m^c4))l/2,  v*  -  [di/di-  d<e!dt  +  <d)/(Acw(z)  +  w/cl  is  the  axial  velocity,  z  — 
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fo  +  /0  (x0»  y0,  £o*  fOdr'  is  the  axial  position  of  the  electron  and  f0  is  the  axial  position  of 

the  electron  relative  to  the  center  of  the  electron  beam  at  t  —  0. 


The  wave  equation  for  the  radiation  field  is  (V2  -  c  2d2/dt2)AR  —  -4 ire  ,  where 

—  \e\ftn  r°°  r°°  r°° 

J  -  — _  dU  J  ^  <b f0  J  "  dy<fi(x 0>  y0)h(t<,) 

-y-lPi8(x  -  Xo)8Cv  -  ^0)8(7  -  z)  (4) 

is  the  current,  (xQ.  y0)  are  the  particle's  transverse  positions  at  t  —  0,  0(xo,  ><>)  IS 
transverse  current  density  profile,  /i(f  0)  is  the  macroscopic  electron  pulse  shape,  nQ  is  the  peak 

I  €  I 

current  density,  and  Pi  —  Aw  is  the  transverse  momentum. 


We  can  rewrite  the  radiation  field  as  A*  —  aR(x,  y,  z ,  t)  exp  [i(a>z/c  —  oif)]e+  +  c.c., 
where  ar  —  AR  exp (/^)  is  the  complex  amplitude  of  the  radiation  field,  and  e±  —  ( ex  ±  iey)/2 
is  a  new  coordinate  system.  The  wave  equation  assumes  the  form 


Ji  +  iU2i-  f+if Lfc,.,„ 

dx2  fiy2  r  I  A*  r  ht  II  * 


c  I  dz  c  Bt 


~  “f  >o)^tfo)  -4h.(2> 

expj-/|j^  (*w(z')  +  <o/c)dz'-  wtjjsfx-  x0)8 (y -  ^0)6 (z -  z). 


In  obtaining  (5),  we  have  used  the  fact  that  Kfitydz2  -  c  2d2/dt2)a,:  I  «  2o»/c|(d/dz  + 


c-'d/dt)fl*|. 


The  solution  for  aR  can  be  separated  into  the  input  radiation  field  a,„,  and  the  excited 
radiation  field,  a„ ,  such  that  aR  —  a,„  +  a„.  The  excited  radiation  field  can  be  shown  to  be  Av 

o„  -  /0  dl'  J  m  d(0  J  ^  dx0  J,  dyof(x0.  >>0,  (0,  x,  y.  z,  r’)8(z  -  c(t- t')  -  z)  (6)  f 

where  !  <L**T**' 


-  -1  x  ,4*(2  -  c(r- r')) 

4ir  c2  (/  -  r')y 


t| m Jr*"1 <*•<*■> + *"> *• -  ' 
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The  integral  in  time  of  Eq.  (6)  can  be  evaluated  by  changing  the  argument  of  the  delta  func- 


fl«  “  Jo  dt  L  d*o  f_„  d*o  f_"  dyof(x0.  y0.  U-  x.  y,  z,  t')  — — 

oz 


*  ,  < 


r*-S0\ix*>***lW  ^  f(x0t  yo,t0,  x.  yt  z,T0) 

Jz~a  °  C-  vz(x0t  yQ,  f0,  r0) 

where  r0(x0,  yo,  t)  is  the  retarded  time  associated  with  the  electron,  which  originated  at 
(*o.  y<*’  ( o)  at  *  “  0-  The  retarded  time  satisfies  the  equation 

fo  +  *c(T0>  +  C(t  ~  T0)  ~  f  +  zc  (/)  (8) 

where  zc(t )  —  v^  (/')<//'  is  the  macroscopic  location  of  the  center  of  the  electron  beam  at 

time  f,  v^fr)  *  &/(kw(zc (r))  +  c o/c)  is  the  macroscopic  velocity  of  the  electron  pulse,  f  is  the 

position  of  the  electron  relative  to  the  center  of  the  electron  beam  at  time  t. 

The  complex  radiation  amplitude  in  Eq.  (8)  can  be  evaluated  if  we  make  the  following 
simplifying  assumptions.  For  experimental  parameters  of  interest,  we  can  assume  that  the 
bunching  mechanism  does  not  alter  the  macroscopic  electron  pulse  shape,  hence,  it  travels 
undistorted  through  the  interaction  region.  We  will  assume  that  the  electron  beam  has  an  axi¬ 
ally  symmetric  Gaussian  profile  in  the  transverse  direction,  i.e.%  0(xo,  .Vo)  " 

exp  l-  (xq  +^02)/r02).  Furthermore,  we  will  assume  that  the  waist  of  the  input  radiation  field  r0 
is  much  larger  than  the  radius  of  the  electron  beam  rby  such  that  i  is  approximately  a  function 
of  and  t  only.  The  excited  radiation  field  takes  the  form 
2  2 

atx(r,  (,  t)  -  -  -j  d(0/i((0) 

,  +  tc(t)  -  c(t- r0))  (1  +  vt/c)y* 

y  cU-T0)-izb  {9) 

Jlk  c(t_  _  — 1  expl-/(i({0.  r0)  +<p(r,  (,  r0))] 
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where  zb  —  rb  o)/2c  is  the  Rayleigh  length  associated  with  the  electron  beam  radius.  Equations 
(3)  and  (9)  describe  self-consistently  a  general,  nonlinear,  2-D,  FEL  amplifier  with  a  macros¬ 
copic  pulse  shape  /f(£0) 


For  the  purpose  of  illustrating  the  finite  length  pulse  effects  in  an  FEL  amplifier  operating 
in  the  low  gain  limit,  i.e.,  |njn|  »  |tfext,  we  take  the  electron  beam  profile  to  be  uniform,  i.e., 
/i(f0)  *  1  for  |f0l  <  LJ2  and  h(£0)  =  0  for  |f0!  >  LJ2,  where  Lb  is  the  length  of  the  elec¬ 
tron  pulse.  We  also  make  the  constant  phase,  resonant  particle  approximation.  In  this  approxi¬ 
mation  all  particles  are  assumed  to  have  the  same  constant  phase,  i  R.  The  electron  beam  in 
this  approximation  consists  of  a  pulse  train  of  macro  particles  separated  in  distance  by  27rv0//a>. 
Furthermore,  we  will  limit  ourselves  at  this  point  to  a  constant  parameter  wiggler  and  consider 
only  an  external  DC  electric  potential.  The  amplitude  and  phase  of  the  total  field  are 


where  Ain 


lr 


AR{rt  t)  ~  A-m  ~  <*o cosjtR  +  /,  sin  i R) 

(r,  f ,  r)  =  —  cxqAhI/,  cos  ijiR  -  lr  sin  4* R] 

-  Re  (/);  /,  =  Im  (/),  /  -  E\~-ql 

I 


(10a) 

(10b) 

—  r2  l  r  . 

— Ei  is  the 

rb 


exponential  integral  function. 


<7/  -  (-&*)( -  f0./)  “  ^by\ 

Qu  *  (-«*)( 2y/*(£  -  Uu)  -fzbV\ 

to.u  *  f  (for  f  <  LJ2)  and  f  *  Lb/2  (for  £  ^  Lj 2)  is  the  upper  limit  of  the  integration,  fo./ 
-  f  — (c—  v*)r  (for  £  -(c  —  v*.)/  >  -LJ 2)  and  £0./  ■*  ~LJ2  (for  f  —  (c  —  v*.)/  <  -LJ 2)  is 
the  lower  limit  of  the  integration,  and  yzR  is  the  resonant  gamma  associated  with  the  axial 
motion. 


A  more  realistic  electron  beam  profile  /t(£0)  *  1  ~(2 (for  f0  ^  LJ 2)  and  Mfo) 
*  0  (for  Ifol  ^  I*/2)  can  also  be  integrated.  The  result  is  not  given  here,  because  the  more 
complicated  expressions  would  obstruct  the  initial  understanding  of  the  physical  process  of  the 
pulse  propagation. 
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As  an  example  of  a  10.6  m  m  FEL  utilizing  a  C02  laser  as  an  input  field,  we  choose  an 
electron  beam  of  energy  25  MeV  (y0  *  50),  current  of  /  *  5  A  and  radius  (Gaussian  profile) 
of  rb  —  0.5  mm  and  pulse  length  Lb  *  3  mm.  Such  a  beam  has  a  peak  density  on  axis  of  n0 
—  1.3  x  1011  cm-3  (a)b  —  2.0  x  1010sec_l).  The  constant  parameter  wiggler  has  a  magnitude 
of  Bw  —  5.0  kG  and  wavelength  of  lw  =  2.8  cm  which  gives  Aw  —  2.2  x  103  statvolts.  The 
wiggler  velocity  is  v0i  *  2.6  x  10“ 2  c  which  gives  yi  *  1.35  and  y2  —  37.  The  input  C02 
power  density  is  taken  to  the  Pin  *  4  x  108  W/cm2  which  gives  Ain  —  0.30  statvolts.  Our  illus¬ 
tration  assumes  resonant  macro  particle  approximation  and  an  applied  D.C.  electric  potential 
such  that  sin  i/i  ^  *■  0.6. 

The  schematics  of  the  gain 

G(r,(,  0-  (AR(r,(.  t)  -  Ain)t Ain 

are  shown  in  Figs.  2  and  3.  The  slashed  bars  in  the  (z.  t)  plot  of  Fig.  2  denote  the  locations  of 
the  electron  beams  at  t\  —  1  m/c  and  t  —  2  m/c,  which  c  is  the  speed  of  light.  The  solid  lines 
in  the  (z,  t)  plot  are  the  light  lines.  The  gain  pulse  on  axis  are  plotted  at  times  and  f2.  We 
see  that  the  excited  radiation  pulse  grows  and  spreads  beyond  the  electron  beam  pulse.  The 
transverse  variation  of  the  gain  at  f  —  0  for  various  times  are  plotted  in  Fig.  3.  The  decrease 
of  radiation  field  far  from  the  axis  is  due  to  refraction  toward  the  center  of  the  beam. 

We  have  obtained  a  general  expression  for  the  growth  of  the  2-D,  stimulated  radiation 
pulse  on  an  electron  beam  of  finite  axial  and  transverse  dimensions  in  an  FEL  amplifier.  We 
included  diffraction  as  well  as  refraction.  In  the  axially  symmetric,  low  gain,  resonant  macro 
particle  limit,  we  have  an  analytical  expression  for  the  radiation  gain.  The  formalism  presented 
here  can  be  modified  to  study  the  radiation  build  up  and  "laser  lethargy"  in  the  FEL  oscillator. 
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